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ABSTRACT 

Extensive surveys of galaxy clusters motivate us to assess the likelihood of cluster- 
cluster lensing (CCL), namely, gravitational-lensing of a background cluster by a fore- 
ground cluster. We briefly describe the characteristics of CCLs in optical, X-ray and 
SZ measurements, and calculate their predicted numbers for ACDM parameters and a 
viable range of cluster mass functions and their uncertainties. The predicted number 
of CCLs in the strong-lensing regime varies from several (< 10) to as high as a few 
dozen, depending mainly on whether lensing triaxiality bias is accounted for, through 
the c-M relation. A much larger number is predicted when taking into account also 
CCL in the weak-lensing regime. In addition to few previously suggested CCLs, we 
report a detection of a possible CCL in A383, where background candidate high-z 
structures are magnified, as seen in deep Subaru observations. 

Key words: cosmology: observations; dark matter; galaxies: clusters: individuals: 
Abell 383; galaxies: clusters: general; gravitational lensing 



1 INTRODUCTION 

The mass density in the central regions of galaxy clusters 
typically exceeds the critical value required for lensing, gen- 
erating multiple-images of background objects. This phe- 
nomenon is known as strong-lensing (SL) and the back- 
ground sources are usually very distant field galaxies, lensed 
into magnified and often multiple arcs on the lens plane. Re- 
cent analyses have shown that many sets of multiply-lensed 
images can be uncovered with high-quality space imag- 
ing measurements and improved modelling techniques (e.g., 
Broadhurst et aL]|2005| Tiesenborgs et a l.||2007[ |Limousin| 
et al. 2008 : INewman et al. 2 009: Zitrin et al.|l2009| |Coe 
et al. | |2010| fpeb et al.|T2010l [Richard et al.||2010| |Merten| 
et al.j2011| ). 

With more precise knowledge of the global and large- 
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scale parameters, and extensive ongoing surveys of galaxy 
clusters in several spectral regions, the possibility of a fore- 
ground cluster lensing a background cluster is of practical 
interest. 

An initial estimate of the possibility of observing 



cluster-cluster lensing (CCL) was made by Cooray, Holder 



& Quashnock ( 1999 1, who predicted that a few dozen CCLs 



may be observed over the full sky. Soon thereafter, two such 



lenses were discovered. Blakeslee ( 2001[ ) and Blakcslce et al 
( 2001[ ) found that the nearby supercluster A2T52 (z = 0.043) 
is actually a chance alignment of two clusters: A2152, and 
a more massive background cluster &t z — 0.134 (which was 
then designated A2152-B). The centres of these two clusters 
are separated by 2.4', and some background cluster galaxies 
of the more distant cluster seem magnified and distorted in 



the image-plane of A2152. Athreya et al. ( 2002 \ have shown 



that an excess of distant galaxies in the South- West area of 
MS 1008-1224 is most likely also a weaker lensing effect of a 
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background cluster near the line of sight. There seem to be 
no other explicit cases of CCLs reported to date. 



are relatively small (and further weakened by the Butcher- 
Oemler effect; 



e.g.. 



Bertin & Lombardi (20011 have also investigated the 



Butcher & Oemler 19781, and in such 



properties of a "double lens" configuration, and mainly its 
effect on WL analyses. In this context, the effect of interest 
is the lensing of a background source by two (at least par- 
tially) aligned lenses, where several such configurations were 



suggested or theoretically discussed before ( 


3.g., Crawford, 


Fabian & ReeSjl986,.Seitz & Schneider,1994 


Molinari, Buz- 


zoni & Chincarinil 1996(1 Wang & Ulmer|1997 


|Gavazzi et al. 


2008). 



Major advances in the capability of detecting weak low- 
brightness emission in the optical and X-ray regions, and 
more recently also in SZ mapping of many clusters, together 
with increased precision in the values of the cosmological pa- 
rameters, make the (strong) CCL phenomenon of practical 
interest as a probe of cluster properties. Additionally, the 
statistics of CCLs enhances the use of clusters as probes of 
the evolution of the large scale structure (LSS). The man- 
ifestation of SL in different bands of the electromagnetic 
spectrum secures the identification of CCLs. Contrasting the 
results from a search of CCLs with theoretical predictions 
may yield important new insight, especially on the late evo- 
lution of the LSS. 

We briefly describe the possible observational signatures 
of CCLs in the optical, X-ray and SZ, and carry out a de- 
tailed calculation of the expected numbers of CCLs in sev- 
eral cosmological models using current values of the global 
and cluster parameters. Our updated treatment here (for 
previous calculation see Cooray, Holder fc Quashnock|1999 1 
yields a wide range of values for the predicted numbers of 
CCLs, reflecting modelling and observational uncertainties. 
Additionally, we report a (possible) discovery of another 
moderately-lensed (magnified by 14 ± 3 %; see sQ, high-z 
background cluster at z ~ O.Q+g l behind A383 {z = 0.19), 
~ 2.2' from its centre, as seen in Figure [2] 

The paper is organised as follows: In §2 we discuss the 
observational properties of CCL. In §3 we present our calcu- 
lation of the probability for CCLs and their predicted num- 
bers. The possible detection of a CCL in the field of A383 
is discussed in 54. Our main results are summarized in 55. 



2 OBSERVATIONAL PROPERTIES OF CCLs 

Lensing of a background cluster results in magnified opti- 
cal images of the background cluster galaxies, and in hith- 
erto undetected signatures in the X-ray and microwave re- 
gions. In the image plane, a clear local overdensity of mag- 
nified, distorted, and stretched optical images would gener- 
ally be expected when the galaxies of a background cluster 
are lensed. The higher redshift of the background cluster 
should result in images fainter by the luminosity distance 
ratio (relative to the lensing cluster), but boosted by the 
magnification effect which though preserves surface bright- 
ness, will magnify the total flux (due to the increased area 
occupied by each source in the image-plane). Also, when 
having multi-band imaging, the higher redshift of the back- 
ground cluster will cause the background galaxies to look 
redder relative to the lens red-sequence galaxies (see Fig- 
ure [2|, though this effect might not be prominent in a sim- 
ple RGB colour-composite image when redshift differences 



a case are more likely to be revealed by producing photo- 
metric catalogues which may exhibit a different, secondary 
red-sequence corresponding to the background cluster. 

The steep dependence of the X-ray surface brightness 
on redshift would generally mean that, even though magni- 
fied by the foreground cluster, the background cluster will at 
best look as a faint part of the foreground cluster emission. 
If the background cluster lies further away from the line of 
sight, there might be a traceable signature, as the magnified 
background fiux will be seen far enough from the foreground 
cluster centre, where the foreground flux is lower and thus 
might enable a clear detection, but only if the background 
cluster is sufficiently luminous. We note, however, that in the 
strong regime this may resemble X-ray images of a substruc- 
ture, merger, or related shocks, and thus without additional 
information, even if such a signal is detected it could well 
be miss-interpreted. In addition (as was noted previously by 



Cooray, Holder & Quashnock 19991, X-ray spectra can be 



used to determine the background cluster redshift, particu- 
larly by the measurement of the relatively strong Fe lines. 
Current measurement capabilities (e.g, with Chandra) en- 
able determining the cluster redshift up to 2; ~ 1. Still, it 
might not be feasible to detect a CCL based solely on X-ray 
imaging measurements. 

The SZ effect is the change in the CMB intensity due 
to Compton scattering of CMB photons as they traverse 



intracluster gas (e.g., Rephaeli 1995 Carlstrom, Holder fc 
Reese|2002[ ). The result is a redshift-independent distortion 



of the CMB spectrum, whose thermal component constitutes 
a decrement of CMB spectrum below ~ 218 GHz, and an 
increment of the spectrum above ~ 218 GHz. 

The SZ effect is measured with respect to the unscat- 
tered CMB at the location of the cluster, so lensing affects 
the SZ signal similarly at all observed frequencies. The to- 
tal observed SZ signal is therefore a sum of the intrinsic SZ 
signal from the lensed cluster, which is then magnified by 
the foreground cluster, plus the SZ signal due to scattering 
in the foreground cluster. 

The fact that the SZ effect is independent of redshift 
may help in making the identification of a CCL more fea- 
sible in SZ surveys than in X-ray surveys. However, it may 
still be hard to disentangle the signals of the foreground and 
background structures, as the SZ signal usually stretches out 
to large (projected) distances and may thus cover-up lensed 
features. The net result from such a trade-off will depend on 
the lens and source redshifts, and on the projected distance 
of the lensed feature from the line-of-sight of the lensing 
cluster, so that generally, structures closer to the line-of- 
sight will be more strongly lensed and magnified, but the 
relative flux from the foreground cluster will also be higher. 
As in the X-ray, also here strong CCL may resemble images 
of merger, or related shocks. Obviously, a CCL identiflcation 
can be more secure if lensing features are revealed in both 
X-ray and SZ measurements, or if clearly detected in optical 
imaging measurements. With photometric and/or spectro- 
scopic data, such a detection could yield precise information 
on both clusters. 

In §4 we elaborate further on the lensing signatures in 
the optical. X-ray and SZ regions, as manifested in A383, 
which seems to be moderately lensing a background cluster. 
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3 PREDICTED NUMBERS OF 
CLUSTER-CLUSTER LENSES 

In order to assess the probability for CCLs we integrate the 
mass function over two cluster populations, namely those 
of the lenses and sources. We do so for both the Press fcl 



Schechter (1974 1 mass function (hereafter PS), and that of 
1 Sheth & Tormen|1999| hereafter ST). For each cluster of the 



lens population we integrate the mass function of sources ly- 
ing behind the lens, and included within the Einstein radius, 
as (properly) determined by the source and lens redshifts. 
The volume element over which the source integration is car- 
ried out is computed in terms of the solid angle defined by 
the Einstein ring in the plane of the lens, projected onto the 
source redshift, by means of the ratio of the squared lens- 
source angular diameter distances. Integrating over the lens 
population then yields the desired number of CCLs. 

Note that we include in this estimate only sources that 
are fully enclosed within the respective Einstein radius, ig- 
noring a partial alignment of the source within this radius. 
In this regard our estimate constitutes a lower limit on the 
number of CCLs. 



3.1 Method 

The equation governing the relation between the concentra- 
tion parameter and the Einstein radius assuming an NFW 
profile ( Broadhurst &: Barkana|[2008 1 is. 



V 3Scr / 



ln(l-fc„)-Ci,/(l + c„) 2;2 



(1) 



where Ac and pi are the overdensity at virialisation and 
critical density at redshift z, respectively, 

2 



Dos 



(2) 



47rG DolDls 

is the critical surface density, Dos, Dol, and Dls are the 
observer-source, observer-lens, and lens-source distances, re- 
spectively, and 

1/3 

(3) 



Rv — 



1.69 



1 + 2 



M IStt^ 



is the virial radius. Also, x = 
radius, and 



g{x) = In ■ 



: tan 



Mpc - h , 
where Re is the Einstein 



(4) 



■■ tanh 




The concentration parameter Cv scales with mass and 
redshift according to the following relation: 

c^=A(M/M,f{l + zf, (5) 

where the parameters A, B, C and Af* are taken from various 
c-M relations as we discuss below. 

The solution of Eq. |[T]) provides the Einstein radius as 
a function of c^, from which the angular Einstein radius 
Be = Re / Da{zi), the ratio between the (physical) Einstein 
radius and the angular diameter distance to the lens, can be 
readily determined. 

Having solved for the angular Einstein radius at the 
lens, we can now estimate the number of source clusters that 




z L 
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Figure 1. The lens-source configuration. The grey area repre- 
sents the section of the lens plane which is included within the 
Einstein radius, as calculated from the lens and source redshifts. 
Integration of the source population is performed along the cone 
extending to the right of the lens, with an angular cross sec- 
tion corresponding to the projected Einstein radius at the source 
redshift, zg. The corresponding source-plane volume element is 
denoted in purple. Note that by virtue of the low angular scales 
involved, one can safely use the flat sky approximation to calcu- 
late the volume element at the source redshift. 



would undergo lensing, i.e. the ones lying behind the lens 
and included within the angular area subtended by the Ein- 
stein radius. For this purpose we set the lens mass and red- 
shift, and integrate the mass function over the relevant mass 
range, a volume element defined by the source-cluster red- 
shift, zl < zs < oo, and the angular diameter element spec- 
ified by the Einstein radius of the lens-source system, pro- 
jected onto the source redshift by means of the squared lens- 
source angular diameter distance ratio, {dAj;^ /dAg )^ (Fig. 1). 
This result provides the number of lensed sources behind a 
lens lying at redshift zl, and having a mass m_t. The total 
number of CCL occurrences is likewise estimated by inte- 
grating the mass function over the mass-redshift space of 
the lens: 



3.2 Results 

As is obvious, results depend significantly on the mass func- 
tion, and quite strongly on the c-M relation. This relation is 
only roughly estimated from numerical simulations of clus- 
ters; therefore, it introduces a large uncertainty in the pre- 
dicted numbers of CCLs. We have used the PS and ST mass 
functions and various c-M relations, each specified in terms 
of a different set of A,B,C and A/, parameters ( in eq. [5|. 
We first use the notation and parameters given i n|Komatsu| 
fc SeljakI (i2002J, ba sed on the work of |Seljak| ( |2000| and 



Bullock et al. (20011. In this notation = 5.2 x IO^^A^g 



(as calculated by us according to WMAP7 parameters) is 
the solution to a{Ai) = Sc, where o-{M) is the present day 
rms mass fluctuations, and Sc the threshold overdensity for 
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spherical collapse at 2 = 0, with A = 10, B = —0.2, and 
C = -1. 

In a ACDM cosmological model with 
{nm,^A,n,h,aa) = (0.266,0.734,0.963,0.71,0.801) as 
taken from WMAP7 results, the above parameters for the 
c-M relation, and clusters in the mass interval, 1 x 10^'' Af© 
- 1 X 10^*^1/0 our calculations yield ~ 0.03 CCLs with 
a PS mass function, and ~ 0.1 CCLs with a ST mass 
function. When taking into account also background groups 
of galajcies down to 5 x IO^^Mq, we obtain ~ 0.1 CCLs 
with a PS mass function, and ~ 0.5 CCLs with the ST 
mass function. 

We repeated the calculation with the c-M relation 



given by Duffy et al. (20081 for their full sample, in which 
M. = 2 xTF^/T^M©, A = 7.85, B = -0.081, and 
C = -0.71. With these values, for a ACDM model with 
(n,„,nA,n,/i,cr8) = (0.258,0.742,0.963,0.719,0.796) taken 
from WMAP5 results (as used in Duffy et al. 2008), and 
mass limits of 1 X 1O"M0 - 1 x our calculations 

yield ~ 1 CCLs with a PS mass function, and ~ 2 CCLs 
with a ST mass function. When taking into account also 
background groups of galaxies down to 5 x 10^^ Mq, our cal- 
culations yield ~ 3 CCLs with a PS mass function, and ~ 8 
with a ST mass function. We note that the c-M relation 



presented in Bullock et al. ( 2001 1 yields similar results 



In order to take into account the lensing projection bias 
of triaxial cluster morphology, the calculation was repeated 
with the c-M relation of Hennawi et al. (20071, in which 
M. = 1.3x lO"/i"^M0, A = 12.3, B = -0.13, and C = -1. 

For these values, in ACDM model with 
{nm,0.A,n,h,a8) = (0.258,0.742,0.963,0.719,0.796) 
taken from WMAP5 results, and mass hmits of 1 X 1O"A/0 
- 1 X lO^'^A^, our calculations yield ~ 9 CCLs with a PS 
mass function, and ~ 17 CCLs with a ST mass function. 
When taking into account also background groups of 
galaxies down to 5 x lO^^Af0, we predict ~ 38 CCLs with a 
PS mass function, and ~ 68 CCLs with a ST mass function. 
Thus, taking into account the lensing projection bias boosts 
CCL numbers by about an order of magnitude. 

We compare these results to the observed c-M relation 
from a small sample of 10 clusters derived by |Qguri et al.| 
( 2009[ ), in which Af . = 1 x lO^^Af©, A = 12.4, B = -0.081, 
and C = —1. This relation yields concentrations higher than 
predicted by ACDM simulations, and even higher than those 
derived observationally by previous work (e.g., [Comerford ^ 
Natarajan 20071, and are likely to be extreme results that 



are perhaps less relevant for our purposes. 

Assuming these values, in ACDM with 
{nm,^A,n,h,a8) = (0.258,0.742,0.963,0.719,0.796) 
taken from WMAP5 results, and mass hmits of 1 X lO"A/0 
- 1 X 10^'' A/0, our calculations yield ~ 165 CCLs with 
a PS mass function, and ~ 260 CCLs with a ST mass 
function. When taking into account also background groups 
of galaxies down to 5 x 10^^ Mq, our calculations yield 
~ 660 CCLs with a PS mass function, and ~ 940 CCLs 
with a ST mass function. Thus, according to the observed 
relation (which is known to produce higher concentration 
than ACDM simulations) the predicted total numbers of 
CCLs are quite large. 

It should be noted that taking into account the weak 
lensing regime, in which the background cluster does not 
have to be within the Einstein radius as projected onto the 



source plane, but can be further out up to several Einstein 
radii, the likelihood of a CCL increases significantly. Specif- 
ically, still assuming the flat sky approximation, the likeli- 
hood of a CCL increases simply as the square of the ratio 
of the projected distance of the background cluster and the 
Einstein radius of the foreground cluster. 

Moreover, extending the calculations to the la ranges of 
the ACDM parameters broadens the ranges of our predicted 
numbers of CCLs by up to a factor of ~ 2. Finally, the use of 
other mass functions - such as |Jenkins et al.| ( |2001[ ) or |Tinker| 



et al. (20081 - can introduce another ~20% variation. 



4 A383: A NEW CCL? 

In the course of this work we inspected (among other clus- 
ters) deep archival multiband images of A383 obtained with 



the SuprimeCam on the Subaru telescope (Miyazaki et al 



2002[ ) in 2002,2005,2007,2008 and new data collected in 2 010 
dedicated to the CLASH sample (Post man et al.||2011 



see 
Ihr 



also below), with total integration times of at least 
and up to ~ 4hr, for each of the B,V,Rc,Ic,i',z' bands. Our 
analysis indicates that this is possibly a new example of 
a CCL, with the lensed system being either a background 
cluster or a group of galaxies. Standard image reduction was 
performed with mscred ta sk in IRAFp] while coadded im- 
ages were created following Nonino et al. ( 2009 1. Zero-points 



were estimated from standard stars observations. 

In these wide-field (~ 30' x 30') Subaru images of the 
field of A383, the multiwavelength coverage uncovers sev- 
eral higher-2 structures. Among these, two large structures 
are seen; one is located ~ 5.2' East and ~ 2.2' North of 
A383, around RA=02:48:24.96 DEC=-03:29:31.8, and the 
second ~ 13' East and ~ 2.8' North of A383, around 
RA=02:48:56.03 DEC=-03:29:06.6, and extends Northwards 
towards a third (possibly different) substructure of similar 
colours (and redshift), but these structures are too far from 
the centre of A383 to be relevant for this work. Our BPZ 
( |Bem'tez||2000| [Benftez et aLlpOOi) |Coe et al.| '2006) pho- 
tometric catalogue, based on the 6 Subaru imaging bands 
mentioned above, suggests redshifts of z ~ 0.3 and z ~ 0.7 
for these structures, respectively, which further reduces their 
expected magnification by A383 to ~ 1. 

The more interesting case we consider here is a clear 
redder sequence of galaxies, ~ 2.5' North-East of the cen- 
tre of A383, around RA=02:48:09.57 DEC=-03:29:41.6 (see 
Figure [2|. The background structure is clearly seen as a red 
overdensity of galaxies very faint in the B-band (see Fig- 
ure |2]for an RGB colour image), implying a higher redshift 
than A383, further confirmed by a photometric redshift of 
z ~ 0.9lg i for all 40 member galaxies. The average pho- 
tometric redshift for these galaxies is 2 = 0.96, though we 
adopt the photo-z of the brightest member, z ~ 0.9, as the 
redshift of this structure. 

We have downloaded Chandra X-ray images and ob- 
tained Bolocam 140 GHz SZ images of A383 to examine 



1 |Valdes| l (T998l l. IRAF is distributed by the National Optical As- 
tronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Figure 2. The central field of A383 with the high-z background cluster marked on the image. Note the much redder colours of the 
galaxies in the background cluster with respect to A383. According to our WL analysis, A383 magnifies the background cluster by 
14 it 3% (see §4 for more details). 



whether this background cluster is seen also in these spec- 
tral regions. As may be expected (see §2), no excess X-ray 
brightness is seen at the location of the cluster, possibly due 
to a low mass (and temperature) cluster or group. Neither 
is the lensed cluster detected in the processed SZ map. The 
optical, SZ, and X-ray signals of the lensed cluster, are used 
below to place limits on its mass, along with a WL analysis 
mass-estimation. 

A383 is one of 25 clusters covered in the CLASH multi- 
cycle treasury program. The CLASH program has been 
awarded 524 orbits of HST time to conduct a multi-cycle 
program that will couple the gravitational-lensing power of 
25 massive intermediate redshift clusters with HST's newly 
enhanced panchromatic imaging capabilities (WFC3 and the 
restored ACS), in order to test structure formation models 
with unprecedented precision. More details can be found in 
[Postman et al.| ( |2011[ ). 

To further base the results of this work, and in the 



framework of the CLASH collaboration, we have constructed 
strong and weak lensing (WL) models for this field (respec- 
tively, Zitrin et al.|2011b Umetsu et al. in preparation). In- 
terestingly, the different background structures are seen in 
the wide-field WL analysis, which we use to derive a mag- 
nification of 14 ± 3 % at the location of the lensed cluster 
discussed here. We note that this cluster is ~ 10 Einstein 
radii away from the centre of A383, therefore the SL model 
cannot be used for the magnification estimate, but only to 
check consistency with our WL analysis of A383 in the re- 
gion of overlap (around 0.7-1 arcminute in radius; see Figure 
5 in |Zitrin et al.|2011b I. 

It should be noted, that these higher-^ structures are 
deduced in both ID and 2D WL analyses, even though dif- 
ferent approaches are employed in these. In the ID analysis 
(e.g., Umetsu et al.|[2"011 1, such background structures are 
seen in the shape of "dips" in the tangential shear profile, 
and translate to local mass excess in the lens convergence 
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Figure 3. Smoothed X-ray image of A383. Ttie location of ttie z~ 
0.9 background cluster is indicated by the red line. No noticeable 
excess is seen in that location, as can be expected (see §2). 



(k) profile (see Okabe et al. 2010 see also Zitrin et al.|201l"b 



Figure 5). In the 2D WL analysis presented here, these will 
simply be seen as 2D structures in the mass distribution, as 
can be seen in Figure [4] 

Note also that the detection here is based on the optical 
data, meaning that the lensed cluster is clearly seen, and not 
by other means which later need further verification. 

In order to calculate, or put limits on the background 
cluster mass, we first examine the total light by this cluster. 
The photometry of the 40 detected z ~ 0.9 cluster members 
add up to a total B-band luminosity of 0.84 



+0.36 



X lO'^L 







(using the LRG templates from Benitez et al. 2009; errors 
include conversion from different bands), in roughly the ~ 2 
square-arcminutes occupied by this structure in the opti- 
cal. After correcting for the 14 ± 3 % magnification effect, 
this yields a total source luminosity of O.TA'^q'^^ x W^^Lq. 



Adopting a typical M/Lb ratio (e.g., Zitrin et al. 2011al, 



of 250±50 M/Lb, this translates into a total mass of 
1.84tQ 5g X for the lensed cluster, which constitutes 

a lower limit on its mass, since we did not account for mem- 
ber galaxies lying below the detection threshold (~26.5 AB 
mag in the i?c-band; 5a). 

Next, we perform a 2D gravitational shear analysis of 
Subaru multiband {BV Rcici' z') data (see Zitrin et al. 2011b 
on A383) to constrain the mass distribution of the z = 0.9 
cluster. The 2D-WL modelling procedure used in this work 



is similar to that of Watanabe et al. ( 2011 1 and Okabe et al 



(20111. More details will be presented in our forth-coming 
paper (Umetsu et al., in preparation). From a photometric- 
redshift (2ph) selected sample of background galaxies {zph > 
0.9, 21 < z' < 26 AB mag), we construct a spin-2 reduced- 
shear field on a regular grid of 0.5' x 0.5' independent pixels, 
covering a 20' x 20' region centred on A383. We model the 
projected mass distribution around the 2 = 0.9 background 
cluster as the sum of two NFW halos, responsible for A383 
and the z = 0.9 cluster, each parameterized with the halo 
virial mass (Mvir), concentration (c„), and centroid position 
(xc,j/c). We use the Markov Chain Monte Carlo technique 
with Metropolis-Hastings sampling to constrain the mass 
model with 8 parameters. 

From a simultaneous 2-component fitting to the 2D 



shear data, we find Mvir = l.5lto.tl x lO^M© (68.3% con- 
fidence limits) for the 2: = 0.9 cluster, where other model 
parameters as well as the source redshift uncertainty are 
marginalized over. With this model, the total projected 
mass within a (typical virial) radius of i? = 1 Mpc is 
M2D — 1.9 X IO^^Mq. This result is in good agreement with 
our optically-based estimation (although the latter consti- 
tutes a lower limit). 

Additionally, we use SZ measurements to deduce the 
gas mass, from which we estimate the total mass. Although 
the z = 0.9 cluster is not detected in the X-ray or SZ im- 
ages, the SZ data do provide a constraint on the gas mass. 
The Bolocam SZ data were processed using the procedure 
described in Sayers et al. (20111; we discuss the relevant 



aspects of this analysis below. In particular, the data are 
effectively high-pass filtered in a complicated and slightly 
non-linear way in order to subtract noise due to fiuctuations 
in the opacity of the atmosphere (i.e., the transfer function 
of the filtering depends weakly on the cluster profile). We 
fit an elliptical generalized NFW profile ( Nagai, Kravtsov & 
Vikhlinin|2007l|Arnaud et al.|2010[ ) to A383, which provided 
a marginal fit quality (xVDOF = 1223/1117) Q 

Using the transfer function computed for this model, we 
then deconvolved the effects of noise filtering of our data to 
obtain an unbiased SZ image of A383. Since the filtering is 
somewhat non-linear, our deconvolved image of the z — 0.9 
cluster will be slightly biased because the transfer function 
used for the deconvolution was determined from the profile 
of A383. However, even in extreme cases the bias is < 10% 
of the cluster peak-height, which is negligible compared to 
our measurement noise for the z — 0.9 cluster. 

In order to estimate the SZ signal from the z — 0.9 
cluster, we first subtracted our best-fit model of A383 from 
our unbiased SZ image. We then computed the integrated 
projected SZ signal within a 1 Mpc aperture (2.14 ar- 
cmin) centred on the z = 0.9 cluster, and found YiMpc ~ 
0.80 X 10~^^ ster, with a statistical error of ay.stat ~ 
1.06 X 10~^^ ster. Restricting to physically allowed posi- 
tive values of YiMpc, this results in a 95% confidence level 
upper limit of YiMpc < 2.66 x 10~^^ ster. Note that the 
model-subtracted A383 SZ fiux within our 1 Mpc aperture 
is YA38:i,modei = 1.89 X 10"" ster. 

Using three published Y — Mgas scaling relations based 
on projected Y (Sayers et al.| |20TT Plagge et al.||2010 and 
|Bonamente et al.||2008j , we estimate the 95% confidence 
level upper limit for the gas mass of the z = 0.9 cluster 
within a spherical region of radius of 1 Mpc to be 5.1, 6.6, 
and 2.7 X 10^'^ M©, respectively. All three of these scaling re- 
lations were constrained largely (or entirely) using clusters 
with much higher masses, which is the likely cause of the 
large scatter in the derived masses. We adopt the gas mass 



limit derived from the Sayers et al. (20111 scaling relation, 



since it was calibrated using Bolocam data analysed iden- 
tically to our A383 data. We note that the spread in mass 



^ Note that we simultaneously fit, and subtracted, a point source 
from the NVSS catalogue at 02:48:22 .09, -03:34:30.5, whi ch we 
found, had a flux density of 12.8 mjy ( [Condon et al.|1998| . This 
point source is approximately 5 arcmin from the centres of both 
A383. and the z = 0.9 cluster. 
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limits from the three scaling relations provides an estimate 
of the uncertainties in the Y/M scaling. 

By assuming a typical gas mass fraction of 13% (e.g., 
|Umetsu"et al._2009j ), we derive a 95% confidence level upper 
limit on the total mass of the z = 0.9 cluster of Msz < 
3.9 X 10^* Mq, in agreement with our previous estimates 
based on the luminosity and our WL analysis. 

For general completeness, we repeat a similar procedure 
also for the X-ray data, although only a rough upper mass 
limit can be expected. In order to estimate the X-ray surface 
brightness from the z ~ 0.9 cluster, we consider X-ray obser- 
vations retrieved from the Chandra archive and carried out 
with the ACIS-I CCD imaging spectrometer, namely obser- 
vation ID 2320, with a total exposure time of approximately 
20 ks, and use for our analysis the best-fit triaxial model of 
the IC gas of A383 obtained by ( Morandi fc Limousin|20TT |. 
The IC gas model generates a theoretical surface brightness 
map which has been subtracted from the raw brightness im- 
age. We then integrate the differential counts within a 1 Mpc 
aperture centred on the z = 0.9 cluster, translated into an 
X-ray luminosity by assuming a metallicity oi Z = 0.3 solar 
value and the photoelectric absorption fixed to the Calac- 
tic value for the coordinates of the background cluster. This 
has been accomplished via fake Chandra spectra, where the 



emissivity model (MEKAL model, Kaastra 1992 Liedahl 



et al. 1995 1 is folded through response curves (ARF and 
RMF) of the ACIS-I CCD imaging spectrometer in order to 
generate theoretical counts for the surface brightness. The 
luminosity value is obtained by rescaling the fake Chandra 
spectrum in order to reproduce the observed number of dif- 
ferential counts. 

Given that the X-ray luminosity mildly depends also 
on the value of the IC gas temperature of the background 
cluster, which is unknown a priori, we start by assuming 
a guessed value of the temperature and derive the corre- 
sponding luminosity. This luminosity is then translated via 
the L—T relation ( Morandi et al. 2007 1 into a global temper- 



ature, which has been used iteratively in order to estimate 
the X-ray luminosity and then the same global temperature 
for the background cluster. This process is repeated until 
convergence of the temperature is achieved. The luminos- 
ity is finally converted into a total mass via standard X-ray 
scaling relations (e.g., Morandi et al. 20071, and errors on 



the physical parameters were calculated via Montecarlo ran- 
domisation. We derive a 95% confidence level upper limit on 
the total mass of background cluster of Mx < 2.5 x lO^^M©, 
in agreement with our previous estimates. 



5 DISCUSSION AND SUMMARY 

We have calculated the predicted numbers of CCLs in the 
ACDM model for different mass functions and cluster prop- 
erties. According to our rather conservative estimates, only 
few (~ 3) CCLs are predicted over the full sky based on 
WMAP7 parameters, using either a PS or a ST mass func- 
tion, for clusters in the mass range IxlO^^-lxlO^'* Mq (for 
both the lens and the source). The number increases some- 
what to ~ 10 when taking into account also background 
groups of galaxies (down to 5 x 10^'^Mq), and considering 
different mass functions, but rises substantially to ~ a few 
dozen when taking into account possible lensing triaxiality 




Figure 4. Contours of the dimGnsionless surface mass density 
K smoothed with a Gaussian of FWHM 1.7', superposed on the 
BRz' pseudocolour image of A383. The lowest contour and the 
contour interval are Ak = 0.03, corresponding to the Icr recon- 
struction error. As can be seen, the WL analysis is sufficiently sen- 
sitive to detect a clear extension around the location of the lensed 
background cluster detected in the optical deep RGB colour in 
Figure |2] 



biases, and to hundreds when considering also the weak lens- 
ing regime. 

Two CCLs were claimed a decade ago - A2152, and 
less significantly, MS 1008-1224, where background galax- 
ies are obviously magnified or create a local over-density 
in the image-plane. In addition, several "double lens" con- 
figurations were suggested or theoretically discussed before 



(e.g., 


Crawford, Fabian & Rees 1986 


Seitz & Schneider 


1994 


Molinari, Buzzoni & Chincarini|1996 Wang & Ulmer 


1997 


Bertin & Lombardi|2001| 


Gavazzi et al.||2008 


I. Comparison 



of the number of observed CCLs with theoretical predic- 
tions (see also Cooray, Holder &l Quashnock|1999 1 is clearly 



important and may add significant new insight on the evo- 
lution of the LSS in ACDM. This could be quite useful in 
light of claimed discrepancies, such as larger than predicted 
Einstein radii, and high concentration or disparities in the 
abundance of giant arcs (e.g., [Hennawi et al.||2007 



Broad- 



hurst fc Barkana[[2008l |Broadhurst et al.||2008| JSadeh fc 



Rcphach 2008 ; Puchwoin k Hilbcrt 2009, Horosh ct al. 2010r 
Meneghetti et al.|2010),Sereno, Jetzer fc Lubini,2010),Zitrml 



et al.|2011a i" 

While inspecting lensing measurements of a sample 
of clusters we have noticed a lensed background cluster- 
ing structure behind A383 (2 = 0.19), in deep Subaru 
imaging. Photometric redshifts imply this overdensity is at 
z ~ 0.9^0' 1, and is clearly seen redder in an RGB colour 
image, and very faint in the B band. Our WL modelling of 
A383 implies a magnification of 14 ± 3 % of this background 
cluster. 

Accounting for this correction, the total B-band source 
luminosity is 0.74lQ 'Jy x lO^^L© summed over all 40 mem- 
bers, which is translated into a lower-limit mass estimate of 
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processed image 




2'48'"15' 48" 0= 47"45' 
RA (J2000) 



deconvolved image 




2'48'"15' 48" 0' 47'"45' 
RA (J2000) 



model-subtracted deconvolved image 




2'48"15' 48^0= 47"45' 
RA (J2000) 

Figure 5. Bolocam SZ images of A383. From top to bottom 
tlie tiiumbnails sliow the liigli-pass filtered (processed) image, 
tlie deconvolved image, and tlie model-subtracted deconvolved 
image. The solid white contours in the processed image repre- 
sent S/N= —1,-2,— 3,.., and the dashed white contours repre- 
sent S/N= +1,-1-2,-1-3,... We do not include S/N contours in 
the deconvolved images due to the significant amount of large- 
angular-scale noise. The white plus sign denotes the centre of the 
2 = 0.9 cluster, and the un-fiUed white circle denotes a radius of 
1 Mpc at z = 0.9 centred on the cluster. The solid white circle in 
the lower left represents the effective FWHM of the PSF in these 
beam-smoothed images. As is in the X-ray, the cluster lies below 
the S/N threshold for detection, allowing only a 95% confidence 
upper mass limit of 3.9 X lO" Mg, see ^ 



l-S^lo.se X IO^^Mq, using a typical M/L ratio of M/Lb = 
250 ± 50. We have also analysed SZ and X-ray data, and 
independent WL measurements, to obtain mass estimates 
of Msz < 3.9 X 10^* Mq, Mx < 2.5 x 10"Mq (95% CL 
upper limits) and Mvir.wL ~ 1.511q 94 x 10^*Mq (or a pro- 
jected mass M2d(< 1 Mpc) ~ 1.9 x lO^^M©), respectively. 
These are commensurate with our prior estimate based on 
the luminosity. These are also in agreement with the fact 
that no excess emission is seen in the SZ or X-ray images of 
A383 at the location of the background cluster (see Figure 
[3|, reflecting its probable low mass. 

Deeper images in the different spectral regions would be 
of interest to examine further the possibility that this clus- 
ter is a CCL, along with spectroscopic data for measuring 
the exact background cluster redshift. In addition, future 
detections of CCLs are important, as their overall number 
could probe the LSS parameters, the cluster (and group) 
mass function, and the c-M relation. 
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